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The broken symmetry that develops below 17.5K in the heavy fermion compound URu 2 Si 2 has 
long eluded identification. Here we argue that the recent observation of Ising quasiparticles in 
URu 2 Si 2 results from a spinor hybridization order parameter that breaks double time-reversal sym¬ 
metry by mixing states of integer and half-integer spin. Such “hastatic order” (hasta:[Latin]spear ) 
hybridizes Kramers conduction electrons with Ising, non-Kramers 5 f 2 states of the uranium atoms 
to produce Ising quasiparticles. The development of a spinorial hybridization at 17.5K accounts for 
both the large entropy of condensation and the magnetic anomaly observed in torque magnetometry. 

This paper develops the theory of hastatic order in detail, providing the mathematical development 
of its key concepts. Hastatic order predicts a tiny transverse moment in the conduction sea, a 
collosal Ising anisotropy in the nonlinear susceptibility anomaly and a resonant energy-dependent 


nematicity in the tunneling density of states. 
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I. INTRODUCTION 

The heavy fermion superconductor URu 2 Si 2 exhibits 
a large specific heat anomaly at To =17.5K signalling 
the development of long-range order with an associ¬ 
ated entropy of condensation, J Q T ° dT ^ ^Rln2 per 
mole formula unit 12 . This sizable ordering entropy 
in conjunction with the sharpness of the specific heat 
anomaly suggests underlying itinerant ordering; how¬ 
ever the anisotropic bulk spin susceptibility (cf. Fig. 
1) of URu 2 Si 2 displays Curie-Weiss behavior down to 
T rsj 70 K indicative of local moment behavior. Initially, 
the hidden order was attributed to a spin density wave 
with a tiny c-axis moment®, which was later shown 
to be extrinsic 5 . However, spin ordering in the form 
of antiferromagnetism is observed at pressures exceeding 
0.8GP#EI. Despite thirty years of intense experimental 
effort, no laboratory probe has yet coupled directly to the 
order parameter in URu 2 Si 2 at ambient pressure, though 
there have been a wide variety of theoretical proposals 
for this “hidden order” (HO) problerrP^. We point the 
interested reader to a recent review on URu 2 Si 2 for more 
details.^ 

Expanding on our recent proposaP^®^ of “hastatic or¬ 
der” in URu 2 Si 2 , here we argue that the failure to ob¬ 
serve the nature of its “hidden order” is not due to its 
intrinsic complexity but instead results from a fundamen- 
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tally new kind of broken time-reversal symmetry associ¬ 
ated with an order parameter of spinorial, half-integer 
spin character. Key evidence supporting this conjec¬ 
ture is the observation of quasiparticles with an Ising 
anisotropy characteristic of integer spin f-moments 22 25 . 
Hastatic order accounts for this unusual feature as a con¬ 
sequence of a spinor order parameter that coherently 
hybridizes the integer spin, Ising f-moments with half¬ 
integer spin conduction electrons; the observed quasipar¬ 
ticle and the magnetic anisotropies thus have the same 
origin. 



FIG. 1: (a) Tetragonal structure of URu 2 Si 2 (b) Tem¬ 

perature dependence of magnetic susceptibility in URu 2 Si 2 

aftePl 


A. Experimental Motivation for Hastatic Order 

Above the hidden order transition, URu 2 Si 2 is an inco¬ 
herent heavy fermion metal, with an large, anisotropic, 
linear resistivity 26 ,and a linear specific heat with 7 = 
~ 200mJmol - 1 K -2 . The development of hidden or¬ 
der results in a significant reduction in the specific heat 
to 70 = ^ 60 mJmol - 1 K -2 , corresponding to the the 

loss of about two-thirds of the heavy Fermi surfacd^. At 
T c = 1.5A, the remaining heavy quasiparticles go su¬ 
perconducting. Under a modest pressure of 0.8GPa, the 
hidden order ground-state of UR^SU undergoes a first 
order transition into an Ising antiferromagnet with an 
staggered ordered moment of order 0.4 /i# aligned along 
the c-axis 7 . de Haas-van Alphen (dHvA) shows that 
the quasiparticles in the hidden order phase form small, 
highly coherent heavy electron pockets with an effective 
mass up to 8.5 mp 2 . Remarkably, these small heavy elec¬ 
tron pockets survive across the first order transition into 
the high-pressure antiferromagnetic phase, leading many 
groups to conclude that the (commensurate) ordering 
wavevectors Q = (0,0,1) of the antiferromagnetic and 
the HO phases are the same® 27 

Perhaps the most dramatic feature of these heavy 
electron pockets is the essentially perfect Ising mag¬ 
netic anisotropy in the magnetic g-factors of the itin¬ 
erant heavy f-electrons in the HO state of UR£P^. This 
Ising quasiparticle anisotropy has been determined by 
measuring the Fermi surface magnetization in an angle- 


dependent magnetic field in the HO state; this magneti¬ 
zation is a periodic function of the ratio of the Zeeman 
and the cyclotron energies, where the former is defined 
through an angle-dependent g-factor g(6) 

AE(d) = g(0)n B \B\. (1) 

Interference of Zeeman-split orbits in tilted fields leads 
to spin zeroes in the quantum oscillation measurements 
(cf. Fig. 2) satisfying the condition 

777 * 

g{9n )— = 2 n+l ( 2 ) 

Tn e 

where n is a positive integer and 0 n is the (indexed) angle 
with respect to the c-axis. Sixteen such spin zeroes were 
identified in the HO state of URSl 22 * 23 1 , and the experi¬ 
mentalists found that 


g c 30 


(3) 


where g _ L = g(6 n ~ f) and g c = g(6 n = 0 ), indicating 
that A E(6 n ) depends solely on the c-axis component of 
the applied magnetic field (T> c ), namely that 


9(0n) = g * cos 6, 


(4) 


where g* = 2.6 in contrast to the isotropic g = 2 for free 
electrons. 

In these high-field measurements, the quasiparticle 
anisotropy manifests itself through the appearance of a 
rapid modulation in the amplitude of the dHvA oscilla¬ 
tions generated by the heavy a pockets of UR^SU as 
the magnetic field is tilted from the c-axis into the basal 
plane. The same magnetic anisotropy is also observed in 
the angular dependence of the upper-critical field of the 
superconducting stat e which develops at low tempera¬ 
tures (cf. Fig. 2)J 24 * 25 l Whereas the dHvA measurements 
could in principle belong to a select region of the Fermi 
surface, the upper-critical field, H c2 {6) is sensitive to the 
entire heavy fermion pair condensate, proving crucially 
that the Ising quasiparticle anisotropy pervades the en¬ 
tire Fermi surface of hidden order state. We note that 
while H c2 (0) matches the anisotropy of the g-factor for 
angles near the c-axis, where H c2 is Pauli limited, when 
the field is in-plane, H c2 (0) is larger than expected, likely 
due to orbital contributions. We note that since the Pauli 
susceptibility scales with the square of the g-factor, 
these resolution-limited measurements of suggest that 

x p (0) = x p *cos 2 e ^7 > 900. (5) 

x± 

Such a large anisotropy should be observable in elec¬ 
tron spin resonance measurements that probe the Pauli 
susceptibility directly in contrast to bulk susceptibility 
measurements where Van Vleck contributions are also 
present. 
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FIG. 2: (a) Anisotropy of the measured g-facW 23 plotted (a) 
in polar co-ordinates derived from spin zeros in quantum oscil¬ 
lation measurements and the anisotropy of the upper critical 
field (b) versus sine of the angle out of the basal plane, show¬ 
ing that the data 23 requires a Pauli susceptibility anisotropy 
in excess of a thousand. 

B. Ground-State Configuration of the Uranium Ion 

Bulk susceptibility measurements (see Fig. 1) of 
URu 2 Si 2 reveal that the magnetic U ions have an Ising 
anisotropy, with a magnetic moment that is consistent 
either with a U 3+ (5/ 3 ) configuration or a non-Kramers 
u A + (5/ 2 ) iodnan. A natural explanation for the quasi¬ 
particle Ising anisotropy is that the Ising character of the 
uranium (U) ions has been transferred to the quasipar¬ 
ticles via hybridization, and this is a key element of the 
hastatic proposal. The giant anisotropy in places a 
strong constraint on the energy-splitting A between the 
two Ising states 

2 A 1 

-RdHvA ) 30 

requiring that in a transverse field, F^hvA = 11T the 
U ion is doubly degenerate to within A ^ IK. Fur¬ 
ther support for a very small A comes from the di¬ 
lute limit, U^Thi^R^SU ( x ± .07), where the Curie¬ 
like single-ion behavior crosses over to a critical loga¬ 
rithmic temperature-dependence 9 below lOiF, log T/Tk, 
where Tk ~ 10 K. This physics has been attributed to 
two channel Kondo criticality, again requiring a splitting 
A <C 10 K. Such a degeneracy is of course natural for 
in a Kramers U (5/ 3 ) ion, containing an odd-number 
of f-electrons. However it can also occur in a 5 f 2 “non- 
Kramers doublet” with a two-fold orbital degeneracy pro¬ 
tected by tetragonal symmetr}^^. 

Motivated and constrained by the bulk spin suscepti¬ 
bility and the quantum spin zeroes, we therefore require 
the U ion to be in an Ising doublet in the tetragonal en¬ 
vironment of URS (cf. Fig. 1 ); such a magnetic doublet 
of URu 2 Si 2 has the form 

|r±) = 5> n |±(J z -4 n)), (7) 

n 

where the addition and subtraction of angular momen¬ 
tum in units of 4 h is a consequence of the four-fold sym¬ 


metry of the tetragonal crystal. However, the presence of 
a perfect Ising anisotropy requires an Ising selection rule 

<r±|J±|iy) = o ( 8 ) 

that, in the absence of fine-tuning of the coefficients a n , 
leads to the condition that — (J z + In') ^ (J z +4n)±l, or 
J z 7 ^ 2(n — n') d= requiring J z G Z must be an integer. 
For any half-integer J z , corresponding to a Kramers dou¬ 
blet, the selection rule is absent and the ion will develop 
a generic basal-plane moment. As an aside we note that 
the fine-tuned case will produce an Ising Kramers dou¬ 
blet, but corresponds to the complete absence of tetrago¬ 
nal mixing, highly unlikely in a tetragonal environment. 

Let us apply this argument specifically to the case of 
URu 2 Si 2 . We first suppose that the U ion is in a 5 / 3 
(J = |) configuration, predominantly in a | ± |) state. 
The presence of tetragonal symmetry results in a crystal- 
field ground-state 

l±> = “I ± =F §> -•- c| =F §> (9) 

so that |(—|J + |+)| 2 = 5 b 2 + 6 ac and perfect Ising 
anisotropy is only achieved with fine-tuning of the tetrag¬ 
onal mixing coefficients such that the condition 5 b 2 + 
6 ac = 0 is satisfied. By contrast, for U ion in a 5 / 2 
(J = 4) configuration, its ground-state may be a non- 
Kramers T 5 doublet 

|±) = a|±3)+6|Tl>, ( 10 ) 

where Ising anisotropy exists for arbitrary mixing be¬ 
tween the | =b 3) and | =p 1) states since this tetragonally- 
stabilized state is dipolar in the c-direction and 
quadrupolar in the a-b plane. Because the phase space 
associated with the non-Kramers doublet is significantly 
larger than that for its finely-tuned Kramers counterpart, 
we take the doublet to be the more natural ground- 
state configuration of the U ion in URS. However we 
have proposed a direct benchtop test to distinguish be¬ 
tween these two single-ion U ground-state configurations 
in URu 2 Si 2 using the basal-plane nonlinear susceptibility 
to check this key assumption in the hastatic proposal 33 . 

The combination of the observed Ising anisotropy and 
tetragonal symmetry are crucial towards pointing us to 
the non-Kramers T 5 doublet. By contrast in a hexagonal 
system, like CeAls 3 -^, the six-fold symmetry mixes terms 
that differ by 6h units of angular momentum, so a pure 
doublet |±M) mixes with |±M') states only if M' = M— 
6n ( n G Z). For J < 7/2, the maximum M — M' is | — 
(—|) = 5 meaning there is no valid choice of M and M'; 
thus there are two Ising doublets for the Ce (J = 5/2) 
case: | ±5/2) and | ±3/2). These Kramers doublets can 
undergo a single channel Ising Kondo effect^ 3 hat will 
differ substantially from the two-channel Kondo physics 
associated with a non-Kramers doublet!^. 
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C. Hybridization, Hastatic Order and 
Double-Time Reversal Symmetry 



FIG. 3: Schematic of (a) conventional (scalar) vs (b) spino- 
rial hybridization where the hybridization is a) a crossover 
and b) breaks spin-rotational and time-reversal symmetries 
and thus develops discontinuously as a phase transition. 


The formation of heavy f-bands in heavy fermion sys¬ 
tems involves the formation of f-resonances within the 
conduction sea. Hybridization between these many- 
body resonances and the conduction electrons produces 
charged heavy f-electrons that inherit the magnetic prop¬ 
erties of the local moments. When this process involves a 
Kramers doublet (the usual case), the hybridization can 
develop without any broken symmetry and thus is asso¬ 
ciated with a crossover. At first sight, the most straight¬ 
forward explanation of hidden order is to attribute it to 
the formation of a “heavy density wave” within a pre¬ 
formed heavy electron fluid. Since there is no observed 
magnetic moment or charge density observed in the hid¬ 
den order phase, such a density wave must necessarily 
involve a higher order multipole of the charge or spin 
degrees of freedom and various theories in this category 
have indeed been advanced. In each of these scenarios, 
the heavy electrons develop coherence via cross-over at 
higher temperatures, and the essential hidden order is 
then a multipolar charge or spin density wave. However 
such multipolar order can not account for the emergence 
of heavy Ising quasiparticles. The essential point here, 
is that conventional quasiparticles have half-integer spin, 
lacking the lack the essential Ising protection required by 
experiment. 

Moreover, in URu 2 Si 2 both optical® and 
tunnelling®®® probes suggest that hybridization 
develops abruptly at the HO transition, leading to 
proposal d 14 * 15 * that the hybridization is an order param¬ 
eter. The associated global broken symmetry and phase 
transition is naturally described within the hastatic 
proposal. As we now describe, hybridization with 
a non-Kramers doublet requires the development of 


an order parameter that breaks double time-reversal 
symmetry, a requirement that leads us to conclude that 
the order parameter has a spinorial quality. (See Fig. [3| 
The observation of heavy quasiparticles with Ising 
anisotropy in the tetragonal environment of URS implies 
an underlying hybridization of half-integer spin electrons 
with integer-spin doublets that has important symmetry 
implications for the nature of the hidden order. More 
specifically such hybridization requires a quasiparticle 
mixing terms in the low-energy fixed-point Hamiltonian 
of the form 


H = \ka)V aa (k)(a\ + H.C. ( 11 ) 

where | ka) and | a) refer to the half-integer spin conduc¬ 
tion and the integer-spin doublet states respectively and 
H.C. is the Hermitian conjugate; here k and a are the 
momentum and the spin components respectively. Be¬ 
cause time-reversal, 0 , is an anti-unitary quantum op¬ 
erator, it has no associated eigenvalue. However double 
time-reversal 0 2 , equivalent to a 2i r rotation, is a unitary 
operator whose quantum number is the Kramers index, 
K = ( — 1 ) 2J where J refers to the total angular momen¬ 
tum of the quantum state; K defines the phase factor 
acquired by its wavefunction after two successive time- 
reversals: 


0 2 |V>) = \^) = K\i{,). (12) 

For an integer-spin state |a), K = 1 since |<a 27r ) = + | ce), 
indicating that it is unchanged by a 2tt rotation. By 
contrast, for conduction electrons with half-integer spin 
states, | ka), \kcr 27r ) = — | ka), so that K = — 1. There¬ 
fore, by mixing half-integer and integer spin states, the 
quasiparticle hybridization 1~L does not conserve Kramers 
index. Indeed application of two successive time-reversals 
to 1~L yields 

e 2 (V\ka)(a\) = V 2n \ka 2n )(a 27r \ = -V 27r \ka)(a\. (13) 

Since the microscopic Hamiltonian must be time-reversal 
invariant, it follows that 

V = -V 2 * (14) 

so the hybridization transforms as a half-integer spin 
state and is therefore a spinor. It then follows that this 
spinorial hybridization breaks both single- and double¬ 
time reversal symmetries distinct from conventional mag¬ 
netism where Kramers index is conserved; we call this 
new state of matter “hastatic (Latin: spear) order”. 

Before we proceed to discuss the theory of spinorial 
hybridization and its consequences, let us pause briefly 
to elaborate on the distinction between spinors and vec¬ 
tors, expanding the previous discussion with emphasis 
on time-reversal symmetry properties. Quantum me¬ 
chanically, the non-relativistic evolution of a wavefunc¬ 
tion i/j(x,t) is described by the Schrodinger equation 
= ihj^r so that 

Q r ip(x, t) — i/j*(x, —t ) 


(15) 
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where 0 is the time-reversal operator; this is an anti¬ 
unitary operation. For vector spins, S -7 —S. The spin 
\ wavefunction is a spinor 


-ft 


so that 


and 


©T(x,£) 




(16) 


(17) 


Q 2 *(x,t) 


-4>l (x,t) 


(18) 


Bosons, with integer spins, are vectors with © 2 = +1, 
whereas fermions with half-integer spins have K = — 1. 
Because of its spinorial character, we can think of 
hastatic order as “hybridization with a twist” since this 
is a simple way of visualization its behavior under 2tt 
(inverted) and 47r (restored) rotations. 

Hybridization in heavy fermion compounds is usu¬ 
ally driven by valence fluctuations mixing a ground-state 
Kramers doublet and an excited singlet (cf. Fig. la). In 
this case, the hybridization amplitude is a scalar that de¬ 
velops via a crossover, leading to mobile heavy fermions. 
However valence fluctuations from a 5/ 2 ground-state 
create excited states with an odd number of electrons 
and hence a Kramers degeneracy (cf. Fig. lb). Then 
the quasiparticle hybridization has two components, 
that determine the mixing of the excited Kramers dou¬ 
blet into the ground-state. These two amplitudes form a 
spinor defining the “hastatic” order parameter 


T = 


*T 

*1 


(19) 


Loosely speaking, the hastatic spinor is the square root 
of a multipole 


T ^ yj multipole. 


( 20 ) 


Moreover, the presence of distinct up/down hybridiza¬ 
tion components indicates that T carries a half-integer 
spin quantum number; its development must now break 
double time-reversal and spin rotational invariance via a 
phase transition. 

Under pressure, URa^SU undergoes a first-order phase 
transition from the hidden order (HO) state to an an- 
tiferromagnet (AFM)P. These two states are remark¬ 
ably close in energy and share many key features 27 29 
including common Fermi surface pockets; this motivated 
the recent proposal that despite the first order transition 
separating the two phases, they are linked by “adiabatic 
continuity,”^ corresponding to a notional rotation of the 
HO in internal parameter space 10 30 . In the magnetic 
phase, this spinor points along the c-axis 




T b 


( 21 ) 


corresponding to time-reversed configurations on alter¬ 
nating layers A and B, leading to a staggered Ising mo¬ 
ment. For the HO state, the spinor points in the basal 
plane 




1 / e - ^/ 2 

71 l e^/ 2 




1 /_ e -*4>/2 

71 V e^ /2 


( 22 ) 


where again, T# = OT A . This state is protected from 
developing large moments by the pure Ising character of 
the 5/ 2 ground-state. 


D. Two-Channel Valence Fluctuation Model 

We next present a model that relates hastatic order 
to the particular valence fluctuations in URv^SU, based 
on a two-channel Anderson lattice model. The uranium 
ground state is a 5/ 2 Ising doubletP, which then fluc¬ 
tuates to an excited 5/ 3 or 5 f 1 state via valence fluctu¬ 
ations. The lowest lying excited state is most likely the 
5/ 3 (J = 9/2) state, but for simplicity here we take it to 
be the symmetry equivalent 5 f 1 state, and assume that 
fluctuations to the 5/ 3 are suppressed - in this sense, we 
take an infinite-U two-channel Anderson model. 

As we can write the f 2 Ts doublet, |±) = a|±3)+6|=Fl) 
in terms of combinations of two J = 5/2 f-electrons in 


the three tetragonal orbitals, Tf and T6, 


1 +) 

I-) 


ftAVr+t' 


■?// 


r et^r+t 




' s /r 6 -7r“-|ft^ 

s/rvi/r-ftl^’ 


(23) 


if we assume a 5 f 1 Ty" excited state, we can now read off 
the valence fluctuation matrix elements directly. Valence 
fluctuations occur in two orthogonal conduction electron 
channels, Tj and T6, and we find 

h V fU) = ^4 6 ±(i)lF±)(r 5 ±| 

+ ^4 7T (j)|r+T)(r 5 ±| + H.c.. (24) 

where d= denotes the “up” and “down” states of the 
coupled Kramers and non-Kramers doublets. The field 


J 

C Fa 


U) = Ek 4(k) 


r»t o * R- 


,j creates a conduction 


electron at uranium site j with spin cr, in a Wannier or¬ 
bital with symmetry T E {6, 7}, while Vq and Vj are the 
corresponding hybridization strengths. The full model is 
then written 

H = Y, *±<L<** + £ [HvfU) + Had)] (25) 

kcr j 

while H a (j) = AEJ2± |r 7 ±,j)(r 7 ±, j\ is the atomic 
Hamiltonian. 

To encompass the Hubbard operators in a field-theory 
description, we factorize them as follows 

x aa = lifter) (r 5 a| = txa■ 


rsj 


(26) 
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E. Structure of the Paper 


{*.) <»£> 



Ising Doublet 


X aa —■► (^^Xa 


FIG. 4: Showing conduction electron self-energy E c . Hy¬ 
bridization with spinorial order parameter (^f a ) permits the 
development of a Ising resonance inside the conduction 
sea, represented by the above Feynman diagram. 


Here | T^a) = xjjfi) is the non-Kramers doublet, rep¬ 
resented by the pseudo-fermions xL while 'Fj. are slave 
bosond^H representing the excited / 1 doublet \T^o) = 
Tj. \Q). Hastatic order is realized as then condensation 
of this bosonic spinor 

^IrXcx o)Xon (27) 

generating a hybridization between the conduction elec¬ 
trons and the Ising 5/ 2 state while also breaking double 
time reversal (see Fig. [4). These slave bosons play a 
dual role in capturing the hybridization while also acting 
as Schwinger bosons describing a 5 f 1 magnetic moment 
with a reduced amplitude, 2 — rif. The tensor product 
Q a p = describes the development of composite or¬ 

der between the non-Kramers doublet and the spin den¬ 
sity of conduction electrons. Composite order has been 
considered by several earlier authors in the context of two 
channel Kondo latticed^ 44 in which the valence fluctu¬ 
ations have been integrated out. However, by factorizing 
the composite order in terms of the spinor we are 
able to directly understand the development of coherent 
Ising quasiparticles and the broken double time-reversal. 

There are two general aspects of this condensation 
that deserve special comment. First, the two-channel 
Anderson impurity model is known to possess a non- 
Fermi liquid ground-state with an entanglement entropy 
of \k,B In The development of hastatic order in the 
lattice liberates this zero-point entropy, accounting nat¬ 
urally for the large entropy of condensation. As a slave 
boson, T carries both the charge e of the electrons and 
the local gauge charge Qj = +x]xj of constrained 
valence fluctuations, its condensation gives a mass to 
their relative phase via the Higgs mechanisnf 46 . But 
as a Schwinger boson, the condensation of T this pro¬ 
cess breaks the SU (2) spin symmetry. In this way the 
hastatic boson can be regarded as a magnetic analog of 
the Higgs boson. 


To recap, here we are arguing that the observation of 
an anisotropic conduction fluid in URu 2 Si 2 indicates the 
coherent admixture of spin \ electrons with integer spin 
doublets, leading us to propose that the order parameter 
in URu 2 Si 2 is spinorial hybridization that breaks both 
single- and double-time reversal. In conventional heavy 
fermion materials, hybridization is driven by valence fluc¬ 
tuations between a Kramers doublet and an excited sin¬ 
glet in a single channel. The hybridization carries no 
quantum numbers and develops as a crossover resulting 
in heavy mobile electrons. However if the ground-state is 
a non-Kramers doublet, the Kondo effect occurs via an 
excited state with an odd number of elections that is a 
Kramers doublet. The quasiparticle hybridization then 
carries a global spin quantum number and has two dis¬ 
tinct amplitudes that form a spinor defining the hastatic 
order parameter 

*=(x) • < 2s > 

The onset of hybridization must break spin rotational in¬ 
variance in addition to double time-reversal invariance 
via a phase transition; we note that optical, spectro¬ 
scopic and tunneling probes in URu 2 Si 2 indicate the hy¬ 
bridization occurs abruptly at the hidden order transition 
in contrast to the crossover behavior observed in other 
heavy fermion systems. 

We now describe the structure of this paper. The mi¬ 
croscopic basis of hastatic order is presented using a two- 
channel Anderson lattice model in section II, along with 
the mean field solution. In section III, we develop the 
Landau-Ginzburg theory of hastatic order, including the 
appearance of pressure induced antiferromagnetism and 
the nonlinear susceptibility, while in section IV, we de¬ 
rive and discuss a number of experimental consequences 
of hastatic order, showing the consistency of hastatic or¬ 
der with a number of experiments, including the large 
entropy of condensation and tetragonal symmetry break¬ 
ing observed in torque magnetometry, and then making 
a number of key predictions, including a tiny staggered 
basal plane moment in the conduction electrons. We end 
with discussion and future implications in section V. 


II. LANDAU THEORY: PRESSURE-INDUCED 
ANTIFERROMAGNETISM 

A. Thermodynamics 

Hastatic order captures the key features of the ob¬ 
served pressure-induced first-order phase transition in 
URu 2 Si 2 between the hidden order and an Ising anti¬ 
ferromagnetic (AFM) phases. The most general Landau 
functional for the free energy density of a hastatic state 
with a spinorial order parameter T as a function of pres- 
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sure and temperature is 

m = a(T c -T)|4'| 2 +/3|'F| 4 - 7 ('F t ( x z 'F ) 2 (29) 


where 7 = 5 (P — P c ) is a pressure-tuned anisotropy term 
and 


T = r 


( cos( 6 > / 2 )e 2 ^/ 2 \ 
Vsin(^/ 2 )e-^/ 2 ; 5 


(30) 


where 0 is the disclination of T^aT from the c-axis, and 
|T | 2 = r 2 . Experimentally the Tafm(P) line is almost 
vertical, indicating by the Clausius-Clapeyron relation 
that the there will be negligible change in entropy be¬ 
tween the HO and the AFM states. Indeed these two 
phases share a number of key features, including com¬ 
mon Fermi surface pockets; this has prompted the pro¬ 
posal that they are linked by “adiabatic continuity”, as¬ 
sociated by a notational rotation in the space of internal 
parameters. This is easily accommodated with a spinor 
order parameter; for the AFM phase (P > P c ), there is 
a large staggered Ising f-moment with 


oc 




(31) 


corresponding to time-reversed spin configurations on al¬ 
ternating layers A and B. For the HO state (P < P c ). 
the spinor points in the basal plane 


Va 



/ g — i <£/2 

( e ^/ 2 


’F.b oc 


V 2 ( 


_g-* 0/2 

,* 0/2 


(32) 


where T# = ©T 5 and there is no Ising f-moment, con¬ 
sistent with experiment, but Ising fluctuations do exist. 

According to the above expression for T, (TV^T) = 
r 2 cos 6 so that the Landau functional can be re expressed 
as 


/ = - a (T - T c )r 2 + (/3 - 7 cos 2 <9)r 4 . (33) 



FIG. 5: Global phase diagram predicted by Landau theory. 


If P < P c , then 7 < 0 and the minimum of the free energy 
occurs for 0 = 7 t/2 , corresponding to the hidden order 
state ordered state. By contrast, if P > Pc, then 7 > 0 
and the minimum of the free energy occurs at 6 = { 0 , 7 r}, 
corresponding to the antiferromagnet. The “spin flop” in 
6 at P = P c corresponds to a first order phase transition 
between the hidden order and antiferromagnet (See Fig. 

§• 


B. Soft modes and dynamics 


The adiabatic continuity between the hastatic and an¬ 
tiferromagnetic phases allows for a simple interpretation 
of the soft longitudinal spin fluctuations that have been 
observed to develop in the HO stat d 47 | 48 J , and even to go 
soft, but not critical upon approaching T#o 49 - These 
longitudinal modes can be thought of as incipient Gold- 
stone excitations between the two phased^. Specifically, 
in the HO state, rotations of the hastatic spinor out of the 
basal plane will lead to a gapped Ising collective mode 
like the one observed. Within our Landau theory, we can 
study the evolution of this mode with pressure. 

In order to study the soft modes of the hastatic or¬ 
der, we need to generalize the Landau theory to a time- 
dependent Landau Ginzburg theory for the action, with 
action S = f Ldtd?x , where the Lagrangian 

= /[vtr] + p (|V\I /| 2 - c- 2 |vi>| 2 ) , 

and p is the stiffness. Expanding T around its equilib¬ 
rium value T 0 , we take = 0 for convenience and write 


4-(x,i) = ijr oe ^(*K/ 2 = (l+i/2^2S6(q)e^' s - ujt a y Ao. 

(34) 

This gives rise to a change in 'lb<r'F = x I’Fol 2 + 
89{x)z\^q\ 2 corresponding to a fluctuation in the longi- 
tudinal magnetization. This rotation in T does not affect 
the first two isotropic terms in /[T]. The variation in the 
action is then given by 


5S = pI'Fol 2 y; \0(q)\ 2 (q 2 - 
q ' 

The dispersion is therefore 


w 2 = (cq) 2 


where 


A 2 = 


2 5{P C - 
P 


c 2 p 

-Wol 2 ) 

(35) 

+ A 2 

(36) 

H^oi 2 , 

(37) 


so that even though the phase transition at P = P c is 
first order, the gap for longitudinal spin fluctuations is 


Aoc I'FoIx/Pc —P. 















Since dP c /dT c is finite, close to the transition, 
\fP c — P ~ \JdP c /dT c (T — T c ), and A oc yjT — T c . In¬ 
elastic neutron scattering experiments can measure this 
gap a function of temperature at a fixed pressure where 
there is a finite temperature first order transition. The 
iron-doped compound, URu^^Fe^SU can provide an at¬ 
tractive alternative to hydrostatic pressure, as iron dop¬ 
ing acts as uniform chemical pressure and tunes the hid¬ 
den order state into the antiferromagnetPF. 

However, higher order terms, like A(TV^T) 2 |T | 2 are 
also allowed by symmetry, and will introduce a weak dis¬ 
continuity of order |To| 4 , 

A 2 = 2<5(Pc ^~ p) |fr 0 | 2 + j\^o\ 4 + O(|frol 6 ). (38) 

So there will likely always be a discontinuity in the lon¬ 
gitudinal spin fluctuation mode at the first order phase 
transition, although it will be of a lower order (|T 0 | 4 ) 
than a generic first order transition (| I) * 


III. MICROSCOPIC MODEL: TWO-CHANNEL 
ANDERSON LATTICE 

Hastatic order emerges as a spinorial hybridization 
between a non-Kramers doublet ground state and a 
Kramers doublet excited state. In our picture of 
URu 2 Si 2 , a lattice of 5 / 2 (J = 4) U 4+ ions provide the 
non-Kramers doublet (T 5 ), which are then surrounded 
by a sea of conduction electrons that facilitate valence 
fluctuations between the 5 / 2 non-Kramers doublet and 
a 5 f 1 or 5 / 3 excited Kramers doublet. The two-channel 
Anderson lattice model has three components: 

H = H c + '£[H VF {j) + H a (j)], (39) 

j 

a conduction electron term, i7 c , a valence fluctuation 
term capturing how the conduction electrons hop on and 
off the U site, Hvf and an atomic Hamiltonian capturing 
the different energy levels of the U ion, H a . 


A. The Valence Fluctuation Hamiltonian 

1. The 5 f 1 model 

The 5 f 2 T 5 non-Kramers doublet is given by 

|5/ 2 :r 5 ±> =o|±3)+6|Tl), (40) 

and all energies are measured relative to the energy of 
this isolated doublet. In principle, valence fluctuations 
may either occur to 5 f 1 (J = 5/2) or to 5 / 3 (J = 9/2), 
and in fact the lowest lying excited state is likely to be 
the 5 / 3 state. However, for technical simplicity, we take 
the lowest lying valence fluctuation excitation to be the 
5 f 1 Ty excited state, 

|5/ 2 :r s ±)^|5/ 1 :r+±)+e- (41) 


where 


I5 / 1 :T+±) = 7?| ± 5/2) + <5| T 3/2) (42) 

is the excited Kramers doublet. The following section 
will show how particle-hole symmetry can be used to for¬ 
mulate the equivalent model with fluctuations into a 5 / 3 
Kramers doublet. 

To evaluate the matrix elements for valence fluctua¬ 
tions we need to express the 5 / 2 state in terms of one- 
particle states. The state can be rewritten as a prod¬ 
uct of the one particle J = 5/2 f-orbitals, |5/2,m) = 
//JU), using the Clebsch Gordan decomposition 


I 1) — (^/il/243/2 + ^4±3/24 v2^ !«>. 

I ±3) = /I 5/2 4 1/2 |fi). (43) 

Next we decompose the one-particle f-states in terms of 
the one-particle crystal field eigenstates, |T±) = /p ± |f}); 
writing = /jL(T/3|ra), or more explicitly, 

tti/2 = /r 6 ± 

4=3/2 = s fl +±- v 4 - ± 

4=5/2 = r l4+± + S 4 ¥± ( 44 ) 

the non-Kramer’s doublet can be written in the form 


|r 5 +> 

|r 5 -) 




' 4/r 6 -|-/r+ 




t 

r+t 
flu- 


sfkdl-^ 1 °) 

(45) 


where p = q = 65 — <277, s = — hS^J | — ap. 

Valence fluctuations from the ground state (5f 2 T$) to the 
excited state (5f x T?) are described by a a one-particle 
Anderson model with an on-site hybridization term 


HvfU) = 


E 


r=r 6 ,r±;a 


vrc{ a (j)fra(j) + H.c 


(46) 


where the vr are the hybridization matrix elements in the 
three orthogonal crystal field channels and c^ a (j) creates 
a conduction electron in a Wannier state with symmetry 
Ter on site j. 

Now we need to project this Hamiltonian down into 
the reduced subspace of the ground-state |5 / 2 : ^a) 
and 1 5 f 1 : T^ a) excited state doublets, replacing 


fra(j) -» ]T \ r P') (T+a'\fMj)\T 5 a) (T 5 a\ (47) 

cr' ,Ot=± \ / 


in (46). Using (45), the only non-vanishing matrix ele¬ 
ments between the 5 f 2 and 5 f 1 states are 


<r+T|/ r - T |r 6 ±) = P , 
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( r 7 ± l/r 6 ±l±) = Q- 


(48) 


Matrix elements for the Ty channel identically vanish 
(T+cr\f r+± \cr') = 0, so the third term in (45) does not 
contribute to the projected Hamiltonian. The final pro¬ 
jected model is then written 


H = £ CkcLPka + E l H ^U) + H a {j)\ 


(49) 


kcr 


where c£ creates a conduction electron of momentum k 
spin cr, with energy ek, and 

HvfU) = Vfe4 8 ±0')|rt±)<r 5 ±| 

+ UcL (i)|r+ T )(r 5 ±| + H.c. (50) 

1 7 + 

describes the valence fluctuations between the T 5 doublet 
and the excited r/~ Kramers doublet. Here Vq = vr 6 q and 
Vj = v r -p while 

1 7 

H a (j) = ASEl r 7±X r 7±| (51) 

± 

is the atomic Hamiltonian for the excited 5 f 1 : 
Kramers doublet. Notice that, as we have projected out 
most of the possible U states, we are working in terms of 
the Hubbard operators, |T 5 ±) and \Tj±) to describe the 
allowed U states. 


The form-factor is then 


[$rk ] a/3 = E 

(Rat at ,&nn} 


g—ik • (Ratai — b-nn) ^ g —ik-ajVAT 


To further develop the valence fluctuation term, we 
need to determine the form factors relating the conduc¬ 
tion electron Wannier states in terms of Bloch waves, 

4« = E 4/3 [^rkE e_ * k • R ’. (52) 

k 

For a single site interacting with a plane wave, these are 
given by [$rk] a/3 = 4/9 ( k ), where 

5/2 ^ 

4/3 ( k ) = E *3m-f(k)(3m- -||5/2m)(m|r, /?) 

m=— 5/2 

5/2 n - 

m=— 5/2 

(53) 

where <a,/3 G ±. However, in URu 2 Si 2 , the uranium 
atoms are located on a body centered tetragonal(bct) 
lattice at relative locations, Rat n = (±a/ 2 , ±a/ 2 , ±c/ 2 ), 
and the correct form factor must respective the lat¬ 
tice symmetries. Our model treats the conduction band 
as a single band of s-electrons located at the U sites. 
Moreover, we assume that the f-electrons hybridize via 
electron states at the nearby silicon atoms located at 
a tv tv = (±a/2, ±a/2, ±z) where z = .371c is the height 
of the silicon atom above the U atoirPH 


4 / 3 ( a Aw)= E e * k ' RiVJV 2 /^/ 3 (a/v/v)- (54) 

{Rat at ,&nn} 


Here, the term e _zk * aNN y^p(& NN ) is the amplitude to 
hybridize with the silicon site at site a a/-a/a while the 
prefactor e _2k * CRnn-b-nn) [ s the additional phase fac¬ 
tor for hopping from the silicon site a at at to the U atom, 
Rat at directly above it. This form of the hybridization 
can also be derived using Slater-Kostef 52 methods, un¬ 
der the assumption that the important part of the hy¬ 
bridization potential is symmetric about the axis be¬ 
tween the U and Si atom. Notice that this function 
has the following properties: [^rk+c]^ = [$rk ] a/3 and 
[<Frk+Q] a ^ = — [^rk] a ^- We should note that this model 
of the hybridization is overly simplified, in that the U 
most likely hybridizes with the d-electrons sitting on the 
Ru site. Such a d-f hybridization can be treated in a 
similar fashion, and is the subject of future work. 


2. The 5 f 3 case 

For simplicity we have discussed the two channel An¬ 
derson model involving fluctuations from a 5 f 2 ground 
state to 5 f 1 (J = 5/2). However, the more realistic case 
involves fluctuations to 5/ 3 , whose low energy states have 
J = 9/2, and are split into five Kramers doublets by the 
tetragonal crystal field, 

|Ia±) = c A |±9/2) + d A |±l/2)+e A |T7/2), 
|ri±) = a\ ±5/2) +b\ =F3/2) 

|4±) = —b\ ± 5/2) + a| =F 3/2), (55) 

where A G (1, 2, 3) labels the three T 6 Kramer’s doublets. 
There are two generic situations: either a T 7 doublet is 
lowest in energy, and the valence fluctuations are then 
determined by the overlap, 


|r 7 ±> = aA4|r 5 ±) + /¥kir 5 T> 


(56) 
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or alternatively, a T 6 doublet is lowest in energy, with 
the relevant overlap, 

|r 6 ±) = a^ 7T \r 5 ±) + p4 t \t 5T ), (57) 

where the form-factors are as above. In both cases fluc¬ 
tuations will involve conduction electrons in both Tq and 
r 7 symmetries. When the lowest excited state is a Ty, 
the valence fluctuation Hamiltonian is given by, 

HvfsU) = Wj r 6 T |r 5 ±)(r 7 ±|+H.c) 

+ W]r 7T |r 5 T)(r 7 ±|+H.c.). (58) 

The only difference between the 5/ 3 and the 5 f 1 model 
is that the excited state requires adding a particle, so the 
valence fluctuation term here is the particle-hole conju¬ 
gate of the 5 f 1 case. 

B. Slave particle treatment 

Hubbard operators cannot be directly treated with 
quantum field theory techniques since they do not sat¬ 
isfy Wick’s theorem. We follow the standard approach 
and introduce a slave particle factorization of the Hub¬ 
bard operators that permits a field theoretic treatment, 

|r+a)(r 5 a| = ^ Xa . ( 59 ) 

where 


We may interpret this replacement as a kind of multi¬ 
particle contraction of the many body Hubbard operator 
into a single fermionic bound-state. Once this bound- 
state forms, a symmetry-breaking hybridization develops 
between the conduction electrons and the Ising 5 f 2 state. 

The dual Schwinger/slave character of the boson Tj. 
representing the occupation of 5 f 1 means when this field 
condenses, it not only breaks the local U( 1) gauge sym¬ 
metry, but also the global SU( 2 ) spin symmetry. How¬ 
ever, it breaks the U{ 1) gauge symmetry as a slave bo¬ 
son, which have been well studied in the context of heavy 
fermions. The U{ 1) phase of the local gauge symmetry is 
subject to the Anderson-Higgs mechanism, in which the 
difference between the electromagnetic gauge field and 
the internal U{ 1) gauge field acquires a masP^. It is this 
process that gives the x fermions a physical charge. The 
combination of the global and local symmetry breaking 
processes means that the hastatic order parameter can 
be thought of as a magnetic Higgs boson. 

With this slave particle factorization, we can reformu¬ 
late 

J2H a (j)^AE *Uj)*<rU)- (64) 

3 3,0 r = ± 

The valence fluctuation term at each site takes the form 

h V fU) = ^4 6± (i)^(i)x±(i) + v4_ T (j)4(j)x±(i) 
+ H.c.. 7 (65) 


l r 5«) = xlM) (60) 

is the non-Kramers doublet, represented by the pseudo¬ 
fermion xtti while is a slave boson 41 representing the 
excited / 1 doublet 


|r+<7> = |n) (6i) 


that carries a positive charge and a spin quantum num¬ 
ber. Condensation of the spin-^ boson then gives rise to 
the hastatic order parameter 


^ = 



(62) 


We now rewrite this expression in terms of Bloch waves 
by absorbing the momentum dependent Wannier form- 
factors into the spin-dependent hybridization matrix, in¬ 
troducing the operator 


V(k,j) = V^ Ve B] + 
where the matrix 



*,0 


( 66 ) 


(67) 


contains the hastatic boson. The valence fluctuation 
term then becomes 


This condensation process that we may now replace the 
Hubbard operator X aa by a single bound-state fermion 

HXa -> (i)Xa- (63) 


HvfU) = E c kA-(k,i)Xa(i)e- ik -^ +H.C. (68) 


Putting these results all together, our model for URu 2 Si 2 is given by 


H = ^ e k4 ff ^ + E[4a(^r,B]+y 7 $ rT -Bp i ) Xa ( J >-‘ k - R i + H.c.' 

kcr k 

+ E A ; (E^(i)^(i) + E^(i)x a (i)-1) ■ 


■A E E 

j,<r=± 


(69) 
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The respective terms in this Hamiltonian describe the 
conduction electrons, the hybridization between the ex¬ 
cited Kramers and ground-state non-Kramers doublets, 
the energy A E of the excited Kramers doublets. The 
second line of the Hamiltonian describes the constraint 
n^(j) + n x (j) = 1 associated with the slave boson repre¬ 
sentation. Finally notice that if we compare (69) with the 
parent Anderson model (46), the original f-annihilation 
operators in the T 6 and T 7 - channel have been replaced 
as follows: 


fr ea (j) ->■ fr 7 a(j) -> (-Bt(TiXj)a. (70) 

So while x± is a slave particle representing the non- 
Kramers doublet, these operators represent composite 
fermions in the two hybridization channels, with all the 
quantum numbers of an electron. 


C. Mean Field Theory for Hastatic Order 

1. Mean Field Hamiltonian: a spinorial order parameter 


The central element of the mean field theory is the 
hastatic order parameter, described by a two-component 
spinor. We consider the following configurations 


/s /g dQ * R? F4>)/2\ 

* = <*->= [ e i(Q-K j+ <t >)/2 )■ 


(71) 


where Q = (0,0, ^ L ), corresponding to a hybridization 
that is staggered between planes, with a spinorial order 
parameter that points in the Basal plane, rotated through 
an angle (j) around the c-axis. Eventually, we shall choose 
< j) = 7 r /4 to provide a 45° rotation of the scattering t- 
matrix relative to the x-axis, an orientation that provides 
consistency with the measured Xxy anomaly in the bulk 
susceptibilit}®. 

Next, in ( [ 66 ] ) we make the substitution 

<V(k, j)) = V 6 $r„(&}) + v 7 $ r - <B>i- (72) 

It is convenient to write (Hj) in the form 

(B)) = 1*1 Us, (73) 


where 


Uj 


,i(Q * Rj +</>)/2 


0 —i(Q • R,-+0)/2 


(74) 


The gauge symmetry of the slave particle representa¬ 
tion allows us to redefine the f-electrons, Xj = ^jXj to 
absorb the spatial dependence of (Hj) into the redefined 
f-electrons ([70]), so that 

BjXj = I^IXi (75) 


and 


B]a lXj = mUjViUfrxj = mn-t)e iQ ' R ’Xj, (76) 

where n = -E (x + y) is the unit-vector representing the 
orientation of the hastatic spinor. The commensurate na¬ 
ture of the wavevector is important, as here we have used 
the fact that = (— \) z o/ c is real. In this gauge, 

the r6 hybridization is uniform while the Tj hybridiza¬ 
tion is staggered. In preparation for our transition to 
momentum space, we write 

V 6 (k) = |*|V6$ 6 (k), (77) 

and 

V 7 (k) = |*|y 7 $ 7 (k)(n-«7), (78) 

where we introduce the short-hand $ 6 (k) = $r 6 (k) and 
$ 7 (k) = $ r _(k). Notice that V 6 (k + Q) = -V 6 (k and 
V 7 (k + Q) = — V 7 (k) both change sign when shifted by 
Q. In the slave formulation, the atomic Hamiltonian is 
H a (j) = ^ ja = AE\^. 

While we treat the hybridization between the conduc¬ 
tion electrons and the f-moments very carefully, we take 
a simple model of the conduction electron hopping, treat¬ 
ing them as s-wave electrons located at the U sites, hop¬ 
ping on a bet lattice with dispersion 


k x a k v a k z c 

6 k = —ot cos - cos —— cos- 11 . 

2 2 2 


(79) 


We do, however, want to capture the essential charac¬ 
teristics of the URu 2 Si 2 bandstructure - namely nesting 
between an electron Fermi surface about the zone center 
and a hole Fermi surface at qR In order to favor a stag¬ 
gered hybridization, and to match up with ARPES ex¬ 
periments suggesting a heavy f-band 5 ^, we take the hole 
Fermi surface to be generated from a weakly dispersing 
X band. This f-electron hopping will be naturally gener¬ 
ated by hybridization fluctuations above Tho > effectively 
where (B^ B) / 0 while (B) =0. A large N expansion 
of this problem would capture these fluctuation effects, 
but is overly complicated for this problem so we put this 
dispersion in by hand, e/k = —Stf cos ^ cos ^ cos 


is a diagonal unitary matrix. 
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So to summarize, our mean-field Hamiltonian is, 

H = Y, ekC 3-Ck. + J> /k + \)xUx*v + K [(A£7 + A )|*| 2 - l] 


[ C k(T C^6 era (k)Xka + V7 era (k)Xk+Qa) 


H.c 


(80) 


where we have dropped the tilde’s on the Xk and A f s is the number of sites in the lattice. We rewrite this Hamiltonian 
in matrix form 


k) 


H = 


T, ( C k> C k+Q’Xk) Xk+Ci) 


/ c k 0 V 6 (k) 

0 e k+Q -V 7 (k) -V 6 (k) 

Vg(k) — V 7 (k) A k 0 

Vd (k) -V 6 f (k) 0 A k+Q J 



K [(A£7 + A)|tt| 2 -A] . (81) 


where we have suppressed spin indices, made the assumption that the Lagrange multiplier A j = A is uniform, 
equivalent to enforcing the constraint on average, introduced Ak = A + e/k, and used the simplification that Q is half 
a reciprocal lattice vector, making V(k + Q) = — V(k), as shown above. When we diagonalize this Hamiltonian, we 
obtain a a set of four doubly degenerate bands, E^. These eigenvalues can be obtained analytically in the special 
case where the band-structure has particle-hole symmetry, but more generally they must be obtained numerically. 


r 


The corresponding mean field free energy is then 


2. Mean field parameters 


3-1 


F[bU] = -^-5>g[l + e 


-PE krl 1 


k,?7 


+U S [ 2 (A£ + A)|lf-A] , (82) 


where (3 = (kriT )~ 1 . In the work presented here, the 
mean field parameters, |T| and A are obtained by numer¬ 
ically finding a stationary point that minimizes F with 
respect to |T| and maximizes it with respect to A. 



20 


25 


30 


35 


FIG. 6: Band structure of the hastatic order is shown in solid 
blue, while the bare conduction (red) and f (green) bands are 
dashed. The parameters used for this calculation are given in 
section IIII C 21 


In order to plot the band structure and calculate semi- 
realistic experimental quantities, we have chosen the fol¬ 
lowing parameters: the internal hastatic angle, <f> = 7 t/4 
is chosen to reproduce Xxy 7 ^ 0 type tetragonal sym¬ 
metry breaking; t = 12.5 meV is taken to match the 
magnitude of Xxy from the torque magnetometry datsP^ 
li/t = —.075 gives the slight particle-hole asymmetry es¬ 
sential to reproduce the flattening of Xxy at low temper¬ 
atures, and has also been adjusted so that /i + A = 0 at 
T = 0 for consistency with the dl/dV calculations (see 
later section); tf/t= —.025 gives a weak f-electron dis¬ 
persion; the crystal field angle £ = .05 is taken to be 
small, as it is in CeRu 2 Si 2 5 -^and NdRu 2 Si 2 ; Vq/V? = 1 
is arbitrary; and finally V 2 /AE = 2 1 is chosen to give 
2 |\I> | 2 « 15% mixed valency. The actual degree of mixed 
valency in URu 2 Si 2 is unknown, with 15% being an up¬ 
per estimate. The band-structure corresponding to these 
parameters is shown in Fig. [ 6 ] 

A plot of |T| and A as a function of temperature is 
shown in Figure [7] B, for these parameters. |T| controls 
the amplitude of the hybridization gap opening at Tho , 
while A controls the location of the hybridization gap 
center in energy. We plot the band structure, both above 
Tho (dashed lines) and at zero temperature (solid lines) 
in the hastatic phase, and the integrated density of states, 
to show how hastatic order opens up a gap above Ep. 

The total density of states is given by 


AO) = 

k?7 

U r d 3 k 

(27r) 3 


= —Im 

7r 




A’k r] id 


(83) 
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where the integral is over the Brillouin zone. The results 
of a numerical calculation of the above density of states 
are shown in Fig. [8]C. The calculation was carried out 
with a discrete summation over momenta, dividing the 
Brillouin zone into 40 3 points and using a small value of 
S to broaden the delta-function into a Lorentzian. 


FIG. 7: Mean field parameters |T| and A as a function of 
temperature. The parameters used for this calculation are 
given in section [III C 2| 



FIG. 8: Density of states in with hastatic order, for the 

region close to the Fermi energy containing the hybridization 
gap. The parameters used for this calculation are given in 
section IIII C 21 


D. Conduction and f-electron Green’s functions 


In order to calculate various moments and susceptibilities, we will require the full 
function, which can be found from the Hamiltonian by integrating out the f-electrons, 


conduction electron Green’s 


[g c (k,iui)} 1 



■ 0 ) 
liO n e k+Q J 


— Vk 



0 

i(jJ n A k _|_Q J 



(84) 


where 


V k 


( v 6k v 7k \ 

\-Vjk -Hik ) ' 


(85) 


Using isospin, t to represent k, k + Q space, we split the conduction electron energy, e k into eok = 5 (e k + e k+ Q), 
eik = \ (ek — e k +q) into the particle-hole symmetric and antisymmetric parts (and similarly with Aok, Ai k ). So now 
we can write the conduction electron Green’s function as: 


[G c (k, iu ;)] 1 = (iw n - e 0 k) - ei k r 3 - E c (k, iuj n ) 


(86) 
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where 


S c (k,zcj n ) = V— 


iio n — Aok + AikT3 


{iuj n Aok)^ 


^ik 


vt 


(87) 


is the self-energy generated by resonant Kondo scattering off the hastatic order. The scattering terms in this quantity 
are determined by two matrices, a diagonal, symmetry-preserving matrix VkV k and a symmetry-breaking VkT 3 V k . 
We decompose these matrices into their channel and spin components as follows: 


= 14 2 c + + (Akc+-<?)Tl, 

V k r 3 v£ = Vic- + Akc— T2, 

where the only non-vanishing components are: 


V 2 

V kc+ 


vi c - = 


Akc+ 

Akc- 




Vk vt 


= 1 E 


Vk-svy, 
VkV^n 
Vkr 3 V^T2 


4 Tr 

= --Tr 
2 


= -Tr 
2 


V 6 kV 6 f k + VrkV^j 
VfikV^ - v 7 k q k 
'(V6kV4 + V7kVt k ). 
V 6 k V’ k - v 7 k v f > k . 


( 88 ) 

(89) 

(90) 

(91) 

(92) 


We note that the non-zero form factor Ak c - which involves a product of hybridization in different channels V 6 k 
and VVk has a d-wave form factor, reflecting the fact that electron scattering off hastatic order parameter breaks 
tetragonal lattice symmetry. The related form factor that describes the c-electron moments, Akc- is a vector in 
spin space which is even parity in momentum space. Under time-reversal, this component changes sign and therefore 
breaks time-reversal symmetry. 

The Green’s function can now be written in the form [Q]~ X = Atq + Br 3 + C • St\ + Dt^ . The eigenvalues are found 
by taking det [t/ c (k, cj)p = 0 , leading to an eighth order polynomial that cannot be generically solved analytically, 
except in the special case of particle-hole symmetry. The four (doubly degenerate) eigenvalues are E^ and are found 
numerically on a grid of k points. Due to the structure of the Green’s function, we can write, 


Q c {iuj n , k) = 


1 


Tinian 


(^ck(iWn)'TO - B ck (ibJ n )T 3 - C ck (iu n ) • (JT\ - £> ck T 2 ) , 


ylck(^C^n) — (^^n ^0k) [(^^n Aok) A^ k J (^^n Aok)bk c _|_ 

B c \^(iuj n ^ = eik [(^n Aok) A lk l AikV k _ 

. —> 

G c k(^^n) — (J'^n Aok) Akc+ 

DckilUJn) = Aik Akc-- 


(93) 


Once is obtained numerically, this structure makes calculations with the conduction electron Green’s function 
fairly straightforward, as we shall illustrate in sections [IVB| and |V C| The f-electron Green’s function can be obtained 
by integrating out the conduction electrons, and is quite similar 


[g f (k,iw)] = (iw n - Aok) - Ai k T 3 - V t — 


lUJrn 


<Tk + <TkT 3 


(iuj n - e 0k ) 2 - e 2 


V. 


(94) 


ik 


The main difference is that the hybridization terms will be of the form TrV k (l, T 3 )VkT a <T 5 , as given below. We note 
that V^f + = V^ c+ = TrV k Vk, and Ak/_ = Ak c - are the same for the c- and f-electrons. There is only one new, 
non-zero term, 


^k/+ = UrV k VkTi<7 = -Ur 


(d k r 7k + r 7 k v 6k ) 


a 


( 95 ) 


that breaks time-reversal symmetry and generally has a d-wave symmetry. 
The f-electron Green’s function is then, 


Q f {iun, k) = 


1 


n„(^ r , E kri ) 


(^/k(iw n )r 0 - B fk (iuj n )T 3 - C fk (iu n ) ■ St\ - D fk T 2 ) , 
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Afk(^n) = (‘ iuj n ~ A 0 k) [fan “ ^Ok ) 2 “ ^ k ] - (icJ n - e 0 k)^ kc+ 

Bfk(itdn) = ^lk [(^7i ^Ok) Gk] ^lk^kc— 

C/k(^n) = (J'^n e Ok)^k/+ 

-^/k £lkA kc _. 


We will examine the k-space structure of these terms, and the related moments, in section V C 


(96) 


E. Particle-Hole Symmetric Case 


For the special case of a particle-hole symmetric dispersion, where (e k+ Q + /i) = — (e k + /i) and e/ k+ q = — e/ k , we 
can solve the Hamiltonian (81) exactly provided A + fi = 0, so that eok = Aok = A = —fi are both dispersionless. In 


fact, the simple c- and f-dispersions we have chosen already satisfy particle-hole symmetry, so the special case A = —fi 
provides a limit where we can obtain analytic results. In this case, eok = Aok = A, and the the determinant of the 
Green’s function can be calculated from (| 86 |. We wish to evaluate the determinant det [cu — H( k)], where H( k) is the 


matrix Hamiltonian given in (81). By integrating out the f-electrons we can factorize the determinant into a product 


of the full-conduction electron determinant and the bare f-electron determinant as follows obtain 

det [uj — ff(k)] = det [—^ c (k, det [cj — Aok — Ai k T 3 ] 

= det [—<? c (k, w) _1 ] ((cj - A 0 k) 2 - Af k ) 2 

where the overall square in the second factor results from the two-fold spin degeneracy and 


det[— Q c (k,uj) ] = det 


(uj - e 0 k ) - ei k r 3 - V k - 


-A, 


Ok 


Ai k r 3 


v; 


: (uj — A 0k ) 2 — A 2 k 

We now impose particle-hole symmetry, setting Aok = A = eok = —/i. For convenience, we redefine 2 


det [uj - iif(k)] = det[-£/ c (k, uj) ^(z 2 - A 2 *,) 2 

z + Ai k r 3 f 

•' lr 


= det z — ei k T 3 — V k 
= V 2 (z)/{z 2 - A ? k ) 2 

where we have multiplied all four rows of the determinant by z 2 ■ 


z 2 - A 2 

z A lk 


O 2 - A? fe ) 2 


A 2 k and have defined 


V 2 (z) = det (z - e k T 3 )( 2 2 - A 2 k ) - V^(z + Ai k T 3 )v£ 


(97) 

(98) 

u> — \. Then 

(99) 

( 100 ) 


Note that since this is a four dimensional determinant, V 2 {z ) is a twelfth order polynomial. Now by employing the 
shorthand V 2 = V kc+ , V 2 = V kc _ , A + = Ak c + and A_ = Ak c -, and substituting VkV k = V 2 + ( A + • <j)t\ and 


V k r 3 V k = V_ + A_T 2 from (89), we obtain 


T> 2 {z) = det z{z 2 - A 2 k - V 2 ) - (z 2 - A 2 k )ei k r 3 - zA + ■ an - AV 2 t 3 - A_Ai k r 2 


( 101 ) 


If we normalize A + = A + h, then the triplet of matrices (71,72,73) = ((n • <t)ti, 72, t 3 ) forms a triplet of anticommuting 
Dirac matrices, ({ 7 *, 7 j} = 2^), which satisfy the charge conjugation symmetry UjiU^ = — 7 ^, (i G [1,3]), where 
U = (J 272 . Since the determinant is unchanged under this transformation, it is unchanged under a reversal 7 \ —>• —7 
If we take the product of the 7 -reversed determinant with itself, the resulting “squared” determinant is then diagonal, 
giving 


V 4 (z) = Det 


- Afk - C 2 ) 2 - (^ 2 - A 2 k ) 2 e 2 k - , 2 A 2 + - A 2 k V _ 4 - A 2 _A : 


lk 


t 4 


And since the argument of the determinant is a diagonal four-dimensional matrix, 

V(z) = z 2 (z 2 - A 2 k - V 2 ) 2 - (* 2 - Afk) 2 e?v - z 2 A 2 s _ - A 2 ,y _ 4 - A 2 _A 2 


t lk 


lk K - 


lk? 


( 102 ) 


(103) 


Now there can be no poles in det[cj — iif(k)] at z = Ai k , so V(z) must have zeros at z = ±Ai k to cancel the denominator 
in det [uj — H( k )] = V 2 (z)/(z 2 — A 2 k ) 2 . We can factorize the determinant as follows: 

V{z) = (z 2 - A 2 k )(z 2 {z 2 - A 2 k - 2V 2 ) - e 2 lk (z 2 - A 2 k ) - 2e lk A lk V 2 ) + z 2 (V* - A 2 + ) - A 2 k (V 4 + A 2 _). (104) 
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Now since D(z = Aik) = 0 is a zero, it follows that 

V*-V* = A 2 + + A 2 _, (105) 

which also follows somewhat nonintuitively from the form-factor definitions, allowing us to rewrite 


V(z) — (z 2 - Aj k ) ^z 2 (z 2 - Aj k - 2V 2 - e 2 k ) + (eikAik - V 2 ) 2 + A 2 ^ 

= (z 2 - Af k ) (z 4 - 2z 2 (f+ -I- -(ef k + Af k )^ + (eikAik - V 2 ) 2 + A 2 ^. (106) 


Therefore, by (99) 



det 

[w - H( k)] = [(w - A ) 4 - 2 a k (w - A ) 2 + 7 k ] 2 

(107) 

where 

a k = Vkc+ + 2 ( e lk + ^lk) 

7k = ( e lkAlk - Vkc-) +A£ C _. 

(108) 

The square in the determinant reflects a two-fold Kramers degeneracy associated with the invariance of the physics 
under a combined translation and time reversal. The energy eigenstates are then determined by the condition 


(w - A ) 4 - 2«k(w - A ) 2 + 7 k = o, 

(109) 

The four bands are then given by, 

Ek v = A ± \Ja-k ± - 7 ^, 

( 110 ) 

where r] = 1 — 4 labels the ++, H—, - 

—b and bands. 



F. Hybridization gaps 


As hast at ic order is, at its heart, a spinorial hybridiza¬ 
tion, the nature of the resulting hybridization gaps is 
essential to understanding the order. In fact, there are 
generically two types of gaps: those that break one or 
more symmetries: translation, spin rotation, crystal or 
time-reversal symmetries, and those that break no sym¬ 
metries. In our mean-field picture, the intra-channel 
gaps, V£ c ± break no symmetries and are proportional 
to the amplitude of the hastatic spinor, (\E^\k), while the 
inter-channel gap, Ak c - breaks all of the above symme¬ 
tries and is proportional to (\E^<r\I/). The role of these 
form-factors as hybridization gap s is especially clear in 
the particle-hole symmetric case, ( |110| ), where their rela¬ 
tive roles may be distinguished. While in the mean-field 
theory, all hybridization gaps will develop at Tho , we 
believe that hastatic order will melt via phase fluctua¬ 
tions, destroying the coherence of the symmetry break¬ 
ing gap, but keeping the symmetry preserving gaps. The 
existence of two types of hybridization gaps that turn 
on at different temperatures can reconcile the number 
of experimen ts tha t find hybridization gaps turning on 
either at T// J 38 * 39 ^ or around the coherence temperature 
T* ~ 50 — 70lj 4Q * 56 l In addition, these hybridization gaps 


will connect different parts of the Fermi surface as the 
intra-channel gaps carry Q = 0, while the inter-channel 
gap connects the folded bands. 

The symmetry-breaking hybridization gap Ak c -, 
shown in Fig. [ 9 ] (a) has an approximate four-fold “d- 
wave” symmetry about a nematic axis n lying in the 
basal plane (Fig. [ 9 ] (b,c)). While only the s quare of 
the inter-channel gap appears in equation (110), it plays 

the same role as the superconducting gap in composite 
J571 


pairing^, and the nodal structure has corresponds to a 
changing phase of the hybridization between c- and f- 
electrons around the Brillouin zone. The modulus of the 
gap, | Akc-1 carries a nematicity, whereby the moments 
of the gap function squared averaged over the Fermi sur¬ 
face, 


( Al c _k a kp ) oc h a hp (111) 

define a secondary nematic director n = (n x ,n y ) of mag¬ 
nitude 


h x ,y OC ^<T x ,y^ (112) 

proportional to the square of the hastatic order param¬ 
eter. The orientation of the nematicity is set by </>, here 
chosen to be 0 = 7r/4. The tetragonal symmetry break¬ 
ing can be tuned by changing the crystal field parameter, 
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£, but cannot be eliminated. However, as the tetragonal 
symmetry breaking is d-wave in nature, it can not cou¬ 
ple linearly to strain and therefore will not lead to a first 
order structural transition. 


IV. COMPARISON TO EXPERIMENT: 
POSTDICTIONS 


A. g-fact or Anisotropy 



FIG. 9: (a) Three dimensional plot showing the symmetry 

breaking component of the hybridization gap, Ak c - for the 
chosen crystal field parameters, £ = .05 and for 0 = 7r/4. 
Orange and blue represent positive and negative values of 
the gap. The gap lobes are oriented in the plane along the 
nematic director h shown in the figure, (b) Top view showing 
the nematic character of the gap function, aligned along the 
nematic director h (c) Side view along the axis of the nematic 
director n showing its nodal d-wave structure. 


The Zeeman energy is determined by the Hamiltonian 


— B ■ M = — Y vT'T'k • B 

k e^BZ 

where ^k = (c k , c k+Q , Xk, Xk+ Q ) T and 



/2/i£>(T 

0 

0 

V 0 


0 0 

2 fiB& 0 

0 gfiiBo z 
0 0 


0 

0 

0 


gf/^BO- 


(113) 


(114) 


where gf is the effective g-factor of the Ising Kramers 
doublet. In a field, the doubly-degenerate energies, |k ga) 
(a = ± 1 ) are split apart so that A E^ v = \E^ — 
= g^(0)B, so the g-factor is given by g^iP) = 

dABwr] • Now we are interested in the Fermi surface 
dB b^ o 

average of the g-factor, given by 


9iP) 


Xlk?7 9\s.rj(0)5(E\^ rj ) 

Sk r] ^(^kr?) 


These quantities were calculated numerically, on a 40 3 
grid, using gf = 2.9 for the effective g-factor of the lo¬ 
cal non-Kramers doublet. The resulting g-factor in the 
z— direction is reduced to g e ff(0 = 0 ) = 2.6 because 
of the admixture with conduction electrons. The delta- 
functions were treated as narrow Lorentzians 5(E) = 
^Im(E — ir?) -1 , where 77 is a small positive number. The 
g-factors at each point in momentum space were com¬ 
puted by introducing a small field SB into the Hamilto¬ 
nian, with the approximation g\^ v (0) = \E^ — E^\/SB. 


B. Anisotropic Linear Susceptibility 


J 


The uniform basal plane conduction electron magnetic susceptibility acquires a tetragonal symmetry breaking 
component in the hastatic phase, given by 


X* 


-(gH B ) 2 T YE Ti K^ C ( k ’ k + Q> iu n )<? v G c { k + Q,k, 


lUJn 


( 115 ) 


iuj n k 

Expanding this in terms of the conduction electron Green’s function, we obtain 

X xy = -( 9 Hb) 2 T W Tr [a x g c (k, iu n )a y g c ( k, f 


, 'ICJrt 


iu> n k 


= -(ffMg) 2 Y 

k?7 


2(B kv - A 0 k )f(E kv ) + (E kv - A ok ) 2 f(E kv ) 


nT7 / ^77 (^k ?7 E] 


k rj' 
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-E 

rj'^rj 


2(fib, - A 0k ) 2 /(£ kT? ) 

(^kr? — ^kr/) Yir]"^y (Ekrj ~ E\zr]" ) 2 


A£ + A^ + 


(116) 


Note that above integral may be positive or negative. 
The functions / and /' are the Fermi function, f{x) = 

(e~ x / T + l) and it’s derivative, f'(x) = df(x)/dx, re¬ 
spectively. The exact nature of the tetragonal symme¬ 
try breaking is determined by the angle of the hastatic 
spinor, </>; when 0 = tt/4, Xxx = Xyy , but Xxy ^ 0, 
but changing 0 can rotate the direction of the tetragonal 
symmetry breaking. 

In our current mean-field approach, all Kondo behav¬ 
ior develops at the hidden order transition, which would 
lead to an entropy of Rlog2 at Tho • Incorporating 
Gaussian fluctuations should suppress the hidden order 
phase transition, Tho , while allowing many of the sig¬ 
natures of heavy fermion physics, including the heavy 
mass and a partial quenching of the spin entropy to de¬ 
velop at a higher crossover scale, T K . The two-channel 
Kondo impurity has a zero-point entropy of - R i og ^pHJGHa 
which should be incorporated into the hastatic phase. 
There is considerable uncertainty in the entropy asso¬ 
ciated with the development of hidden order, S(Tho), 
due to difficulties subtracting the phonon and other non¬ 
electronic contributions, leading to estimates ranging 
from . 15-R log to .SRlog^l If we take a conservative 
estimate of S(Tho ) = .2i?log2, and the normal state 
7 = 180mJ/molK®, S(T K ) = ,2i?log2 + f£* o 7 dT = 
\R log 2 yields TV = 27K. much lower than the coher- 
ence temperature seen in the resistivity. 

V. COMPARISON TO EXPERIMENT: 

PREDICTIONS 

A. Resonant Nematicity in Scanning Probes 

To calculate the tunneling density of states, we assume 
that the differential conductance is proportional to the 
local Green’s function on the surface of the material 

T-(x) oc A(x,eV) (117) 

where 

A(x, uj) = — ImG acr (x, u) — iS) 

7r 


Vv(x) = Xkc(l x - R il)’^a( x 

3 


f'OO 

= E / 0)})e*qil8) 


is the imaginary part of the local electronic Green’s func¬ 
tion. 

To calculate this quantity, we decompose the local elec¬ 
tron field in terms of the low energy fermion modes of 
the system. Typically, in a Kondo system, there are two 



FIG. 10: Schematic illustrating the tunneling and co¬ 

tunneling into three channels: a conduction, a T6 and a T6 
f-electron channel. 


channels - a conduction channel, and a f-electron chan¬ 
nel into which the electron may tunneP^. However, in 
URu 2 Si 2 the presence of a non-Kramer’s doublet now in¬ 
volves two f-tunneling channels - the T 6 and V 7 - channel, 
and tunneling through these three channels can mutu¬ 
ally interfere (see Fig. [l0| ). We therefore decompose the 
electron field Vv(x) in terms of a conduction and two 
f-electron channels, writing, 


■^j)?0aa( X -^d)) ' [ fjF 6 a ] 

\fjr 7 a ) 


( 119 ) 
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FIG. 11: Upper panel: showing the total density of states, 
decomposed according to f, conduction and total spectral 
weight. These curves were calculated for a particle-hole sym¬ 
metric dispersion, setting A = —fi as described in III C2 
and HIE. The red curve shows energy dependent nematic- 
ity. Small panels below show the density of states at three 
different energies, showing the energy dependence of the ne- 
maticity. Reprinted from 19 . 


where 0 c (|x — Ry|) is the wavefunction of the conduction electron centered at site j, while 

Ox-Rtfl) = ^(Ix-R^Jx-Ri) 

<Ox-R;l) = ^ 7 (I X — Rjl)rJo:(x — Rj) 

are the wave functions of the Tq and T 7 - f-orbitals centered at site j. 

Projected into the low energy subspace, following eqs (70) we have fjr 6 a and fj p 7 

Writing B j = bUj , and Xj — ^jXj the expression for the electron field operator becomes 


lM x ) = E (^c(|x-Rj|),</>® a (x-Rj),^>L(x-Rj)) • 


F/cr 


bXjot 

yb(h • 


Xjc 


Next, rewriting the field operators in momentum space, 


c ja = E e ik ’ R;i (c ko - + e iQ ' Rj c k + Q(J ) 

k e^BZ 

Xja = E e * k R ' ( Xka + e * Q Rj Xk+Qa) 
k e^BZ 


( 120 ) 

( 121 ) 


( 122 ) 
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we can decompose the electron field operator as the dot product of two four component vectors 


Vv(x) 


E e H ’-V/(x R,)- 

j^\BZ 


( Cka \ 

^k+Qa I 
Xka I 
\Xk+Qa/ 


(123) 


where 


a j{x) = (<M|x|)6aaX Q ' R ^c(|x|)<W, &<#L( X ) + e iQ '&^L( X )(“ V), e iQ ' R ^0® Q (x) + 6^ a (x)(n• a)) (124) 


We choose a layer where e hQ* R j) = -j-i ? then on this 
layer the local Green’s function is given by 

G(x,w) = E A(x — Rj) • Gji(w) • A f (x — R;) (125) 
where 

Qji M= E ^[(l + ri)a(k, W )]e- <k -^- R ') 

k e\BZ 

is a trace only over the momentum degrees of freedom, so 
Qji is a four by four matrix for each pair of lattice points 
j and Z, where 

A(x) = (</> c (|x|)<5 CTa ,&<^ a (x) + b<j>l a (-x.)(n-a)) . (126) 

The final spectral function is then 


A(x, uj) = —Im Tr 

7r 


E A(x - R j)- Gji{u> - id) • A f (x - Ri) 
o,i 


(127) 

To evaluate this quantity, the summations were limited 
to the four nearest neighbor sites at the corner of a pla- 
quette. The positions x were taken to lie in the plane 
of the U atoms. The wavefunctions 0 6 (|x|) = e - M/ a , 
0 7 (|x|) = e - l x l/ a and 0 c (|x|) = e - l x l/ a were each taken 
to be simple exponentials of characteristic range equal to 
the U — U spacing a. 

The nematicity of the tunneling conductance was then 
calculated numerically from the spatial integral 


/ A(x, eV) sgn(xy)dxdy 

v{eV) = 7 - 77172 - 

(f dxdyA(x, eV) 2 — [J dxdyA(x, eV)] J 

V (128) 
This nematicity is shown in Fig. El Note that the ne¬ 
maticity near the Fermi energy is nearly zero, consistent 
with the absence of quadrupolar moments, but is largest 
at the hybridization gap energy. 


B. Anisotropy of the Nonlinear Susceptibility 
Anomaly 

1. Landau Theory 


tery. It has been understood phenomenologically within 
a Landau theory as a c onseq uence of a large T/ 2 B 2 cou¬ 
pling of unknown origirJ 63 ®, which can now be under¬ 
stood within the hastatic proposal. While the conduc¬ 
tion electrons couple isotropically to an applied field, 
the non-Kramers doublet linearly couples only to the z- 
component of the magnetic field, B z = BcosO , which 
splits the doublet as it begins to suppress the Kondo ef¬ 
fect. When we include the effect of the magnetic field in 
the Landau theory, we obtain 

m = [a(T c -T)-rj z Bl- V± Bl]^ 2 +l3\^\ 4 +^a z ^f, 

(129) 

where the coefficients of the T 2 L> 2 and the terms, 

r] z and 77 j_, will be estimated using a simplified micro¬ 
scopic approach discussed shortly in Section |V B 2[ Min¬ 
imizing this functional with respect to T, we obtain 


/ = -T [a(T c — T)— r] z (B cos 9) 2 - r, ± (B sin 8) 2 } 2 . 

_ (130) 

Following the arguments oF® we can calculate the jump 
in the specific heat A C v and the linear and nonlinear 
susceptibility anomalies and Ays respectively, to 

find 


A Cy 
Tho 
dx l 
dT 


AX3 


2d 


(131) 


- Jp {vz cos 2 9 + rj± sin 2 0) 

(132) 

^ ( r\ z cos 2 0 + r]± sin 2 #) 2 « cos 4 $133) 

/j /j 


where and Ay 3 are the anomalous components of 

the linear and nonlinear susceptibilities that develops at 
Tro • These results show that Ay 3 will exhibit a giant 
Ising anisotropy; we note that these results are compati¬ 
ble with the observed Fermi surface magnetization results 
that indicate that g(0) ~ cos# so that xi ~ cos 2 The 
thermodynamic relation 


AC 

-7j=rX 3 = 12 



(134) 


The origin of the large c-axis nonlinear susceptibility 
anomaly in URu 2 Si 2 63 has been a long-standing mys- 


is maintained for all angles #; the important point here 
is that the anisotropy in Ay 3 is significantly larger than 
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that in A%i and numerical estimates will be discussed 
once we have introduced the microscopic approach to 
hastatic order. 


2. rj z and g± from Microscopies 

To complete this simple Landau theory, we will calcu¬ 
late ff in a simplified model: we will neglect the momen¬ 
tum dependence of both the f-level and the hybridiza¬ 
tion and take the hastatic order to be uniform. None 
of these assumptions qualitatively changes the results. 
The |T| 2 coefficient is calculated from the microscopic 
theory (see the next section) by expanding the action, 
S = -Trlogfl- in where J 0 = 

(iUn-X-gfUfBzCiz)- 1 and Q 0 = (iw n -e k -g/2B • a) -1 


are the bare % and conduction electron Green’s functions 
(remember, x are the fermions representing the non- 
Kramers doublet). V represents the hybridization ma¬ 
trix elements, which are momentum-independent here, 
and proportional to the unit matrix. Note that while the 
conduction electrons are isotropic, the y’s are perfectly 
Ising. The coefficient of |T| 2 is then, 



= V 2 tVV T---t-^-(,135) 

lU) n 6k cr lUJ n X(j 

iuj n k <j 


where e kcr = e k — g/2<r • B and A a = A — gfPfaB z are the 
dispersions in field. Performing the Matsubara sum, we 
obtain 


V J 


/ oo 

deV{e) 

-OO 


tanh — tanh 4^ 


2 T 


2 T 


2(A a - 6 kcr ) 


= pM 


?/. 


D tanh 

de - 

—D 2(A, 


e—g/2cr • B 
2 T 


tanh 


6+| ff-B) 


(136) 


where we approximated the conduction electron density of states as a constant, p within the bandwidth, 2D. Let us 
first calculate , quantizing the field along the x-direction and taking g = 2, 


V±= v 


V =~pv‘ 




D d 2 tanh e — tanh ^ 


_Y) dB\ 2(A — 6 + a B j_) 


(137) 


B± = 0 


As the integrand is a function of e — aB, the integral is straightforward. And as T > A, T, the dominant term will 
be: 


V± = ~pv 


seen 

4T(e - A) 


+ 


tanh — tanh ^ 


(6-A)2 


D 


-D 


pV 2 


(138) 


r\ z will have three contributing terms: one purely from the conduction electrons that is g_\_, one arising from cross 
terms between the conduction and f-electrons, and finally one solely from the f-electrons that dominates the other 
two. We shall focus on this last term, 


Vz = V 




d 2 tanh ^ — tanh A 9 ^ aBz 
dB 2 z 2(A - e - g f p f (jB z ) 



This integral cannot be done analytically at finite temperature, so we take T —»• 0. 


Vz 


-PV 2 


d 2 ^ r° i 

2(A a — e) 


B z =0 


pV 2 

"A2“ 


pV 2 

{D + A) 2 


pV 2 

+2 ' 
1 HO 


as A = Tho at zero temperature. So rj±/Vz 



(139) 


(140) 


Using a conservative value of Tho/D ~ 1/30, we pre¬ 
dict an anisotropy of about 900 in d\i/dT and nearly 
10 6 in A% 3 . However, in a realistic model, there 
will be f-electron contributions to involving fluctu¬ 
ations to excited crystal field states that may reduce the 


anisotropy somewhat. The important point here is that 
the anisotropies will be orders of magnitude larger than 
the single ion anisotropy in xi (approximately 3), and 
furthermore, that they will develop exclusively at the 
hidden order transition. 
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C. Basal-Plane Moment 

Another key aspect of the hastatic picture is the pres¬ 
ence of broken time-reversal symmetry in both the HO 
and AFM phases, manifested by a staggered moment of 
wave vector Q = (0, 0,7r) 

to(Q) = to c (Q)+to/(Q). (141) 

fh contains two parts: a conduction electron component 


According to the Clogston-Anderson compensation 
theorem^, the magnetic polarization of the conduction 
electrons rh c ~ 0 {Tk/D ) is small and set by the the 
same ratio Tk /D that determines the g-factor anisotropy. 
The magnitude of the second / 3 term is set by the over¬ 
all magnitude of T, which in turn is determined by the 
overall amount of mixed valent admixture of 5/ 3 config¬ 
uration into the ground-state. 

Writing out the conduction electron polarization m c in 
detail using Eq ( [93] ) we have 


™c( Q) = ^Y^E^ C k+Qa^ C k/ 3 +H.c) 

k 

= - 9J f^EE Tr [<r<5 c (k, ioj n )Ti] ■ (142) 

icon k 

which involves the off-diagonal component of the conduc¬ 
tion electron Green’s function (see section mm and an 
f-electron component, 


rhf( Q) = mi(Q)z + mjs(Q). 


(143) 


Here 

TO/( Q) = f T t EE Tr [a z G f (\s.,iu n )Ti] (144) 

iuj n k 

is the Ising 5/ 2 contribution, and 

to/s(Q) = (145) 

is the contribution derived from valence fluctuations into 
the 5/ 3 Kramer’s doublet, where T is the staggered com¬ 
ponent of the hastatic order parameter. In the antifer- 
romagnet, the hastatic order parameter at site j is given 

by 


rn c { Q) = -( gjJL B )T ^ Tr [vG c (k, iu n )ri} 

k ,UJ n 

(-^k?7 ^0k) 


= -(flMs)E 


k?7 


Ekv') 


f{Ekn) A kc +- 
(151) 



T(K ) 


FIG. 12: Predicted temperature dependence of the basal 
plane moment. Parameters used for this calculation are given 
in section. 


where 


T. = exp 


-.(Q-R^fj 




(146) 


4/ = (^’) . (147) 

points out of the plane. By contrast, in the hidden order 
phase 


where 


^i-exp <(Q-Rj)< 





( e^/ 2 


(148) 


(149) 


lies in the basal plane, and (j) determines the angle of 
moment from the x-axis in the plane; now the magnetic 
moment lies entirely in the basal plane, determined by 


m(Q) = m c (Q) 


gp^B 




(150) 


Fig. (12) shows the temperature dependence of the 
magnetic moment calculated for a case where D/Tk ~ 
30, for which =0.015/^, which is an upper bound 

for the predicted conduction electron moment. Neutron 
scattering measurements on URu 2 Si 2 have placed bounds 
on the c-axis magnetization of the f-electrons using a mo¬ 
mentum transfer Q in the basal plane. Detection of an 
m_L(0) carried by conduction electrons, with a small scat¬ 
tering form-factor requires high-resolution measurements 
with a c-axis momentum transfer. Recent high resolution 
neutron measurements to detect this small transverse 
moment have not detected a signal, and quote a bound 
on the f-component of the moment \fhf\ <0.001 fis- We 
discuss the implications of this result in section |VIB 


We note that there have been reports from /iSR and 
NMR measurement^^ of very small intrinsic basal 
plane fields in URu 2 Si 2 comparable with this bound. 

We can also examine the quadrupolar moment associ¬ 
ated with the HO state. This is set by the expectation of 
the transverse components of the non-Kramers doublet, 


Q: 


x,y ^ (Xa( a ^,y)a^X/3) 


2 
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= -T^2 Tr [ (7 x,yG f ('k,iUn)T l] . (152) 

iu) n k 


a) 


If we expand this using the f-Green’s function from eq 
(96), we find 


Q., 


(153) 


The f-electron quadrupole moment (153) has a n identical 
form to the conduction electron moment, (151), with e 
A everywhere, and the relevant form-factor is (95) 


t kf+ = —= —2Tr [ (v 6 f k V 7k + V+ k V 6k ) 


(154) 


which has a d-wave form-factor. This means that the 
summation over momentum vanishes, so that the stag¬ 
gered quadrupolar moments Q x , y must vanish. (Indeed, 
there would be no associated lattice distortion, even for 
a uniform hastatic order. ) As a d-wave quadrupole 
is an L = 4-tupole, or “hexadecapole”, this means that 
like Haule and Kotliai^, the hastatic order has stag¬ 
gered ( J x J y + J y J x )(J% — Jy ) hexadecapolar moments. 
However, unlike Haule and Kotliar, where the hexade¬ 
capolar moments are the primary order parameter (and 
thus of order one), here the hexadecapolar moments are 
a secondary effect of the composite hastatic order, and 
like the conduction electron moments, will be of order 
Tk/D. Given how difficult it is to observe large hexade¬ 
capolar moments, the hexadecapolar moments associated 
with hastatic order will almost certainly be unobservably 
small. By contrast, in the antiferromagnetic phase, the 
f-electrons develop a large c-axis magnetic moment. 

We will return to the predicted basal plane moment 
in the HO phase of URu 2 Si 2 in Section VI B. when we 
discuss recent experimental constraints. 


VI. DISCUSSION AND OPEN QUESTIONS 

In summary, the key idea of the hastatic proposal for 
hidden order in URu 2 Si 2 is that observation of heavy 
Ising quasiparticles implies the development of resonant 
scattering between half-integer spin electrons and integer 
spin local moments. It is perhaps useful to contrast the 
various staggered multipolar scenarios for the hidden or¬ 
der with the hastatic one proposed here. In the former, 
mobile f-electrons Bragg diffract off a multipolar density 
wave (see Fig[l3](a)), whereas in the latter, the multi¬ 
pole contains an internal structure, associated with the 
resonant scattering into an integer spin f-state. (see Fig 
[l3] (b)). Hastatic order can thus be loosely regarded as 
the “square root” of a multipole order parameter, 

T ~ i/multipole OP. (155) 




k + Q 


b) __ 

Integer spiiT ^^ k + Q_ 

xp = -y/Multipole = Spinorial OP 



FIG. 13: Schematic contrasting the multipolar and spinorial 
theories of Hidden order, (a) in a multipolar scenario, the 
heavy electrons Bragg diffract off a staggered spin or charge 
multipole (b) in the hastatic scenario, the development of a 
spinor hybridization opens up resonant scattering with a an 
integer spin state of the ion. The multipole is generated as 
a consequence of two spinorial scattering events. In this way, 
the Hastatic spinor order parameter can be loosely regarded 
as the square root of a multipole. 


nematic (see Fig. 9.), with a director ft = (n x ,n y ) of 
magnitude determined by the square of the hastatic or¬ 
der parameter, 


( n x + in y ) oc (156) 

It can also be viewed to result from a symmetry-breaking 
Kondo effect between non-Kramers and Kramers dou¬ 
blets. Hastatic order should be present in any f-electron 
material whose unfilled f-shell contains a geometrically 
stabilized non-Kramers doublet, and we expect its re¬ 
alization in other 5/ uranium and 4/ praseodymium 
compounds. Praseodymium compounds are particularly 
promising tests for hastatic order, as the presence and 
nature of any non-Kramers doublets can be determined 
via inelastic neutron scattering. Any non-Kramers dou¬ 
blet Pr compound must order either magnetically or 
quadrupolarly or form hastatic order - there is no non¬ 
symmetry breaking option, as in Kramers materials. 


A. Broader Implications of Hastatic Order 

At a microscopic level hastatic order demands a new 
kind of particle condensation, one that gives rise to a Lan¬ 
dau order parameter that transforms under half-integer 
spin or double-group representations! 21 ^ Conventionally 
Landau theory in electronic systems is based on the for¬ 
mation and condensation of two-body bound-states. For 
example the development of a magnetic order parameter 
M(x ) is given by the contraction 


In fact, as we have seen the square of the hastatic or- 1 - 1 

der parameter breaks tetragonal symmetry, and is thus a a/ 3 • M(x) 


( 157 ) 
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and s-wave superconductivity is based on the formation 
of spinless bosons 

W(1)V’J.(2)= -F( 1 - 2), (158) 

where F(1 — 2) = — (74/q( 1)^(2)) is the anomalous 
Gor’kov Greens function that breaks the gauge system 
of the underlying system ( see Fig. [M] a). The take- 
home message from conventional two-body condensation 
is that when the two-body bound-state wavefunction car¬ 
ries a quantum number (e.g. charge or spin), a symmetry 
is broken. However under this scheme, all order param¬ 
eters are bosons that carry integer spin. 



^U)^(2) = -F(l-2) 


b) 



IR 


^(l)V>t(2)V>(3) = G(l, 2,3; R)(&(R)) X (R) 


FIG. 14: Schematic Feynman diagrams indicating (a) two- 
body (b) and three-body electronic bound-states where in the 
latter case spin indices have been suppressed for pedagogical 
simplicity. 


all order parameters are necessarily bosonic and bosons 
carry integer spin. However the connection between spin 
and statistics is strictly a relativistic idea that depends 
on the full Poincare invariance of the vacuum. This 
invariance is lost in non-relativistic condensed matter 
systems,where high energy degrees of freedom are inte¬ 
grated out, suggesting the possibility of order parameters 
with half-integer spin that transform under double-group 
representations of the rotation group. Spinor order pa¬ 
rameters involving “internal” quantum numbers are well 
known in the context of two-component Bose-Einstein 
condensates. The Higgs field of electroweak theory is 
also a two-component spinor. However in neither case 
does the spinor transform under the physical rotation 
group. Moreover it is not immediately obvious how such 
bound-states emerge within fermionic systems. 

Hastatic order is a generalization of Landau’s order 
parameter concept to three-body bound-states. This is 
natural in heavy fermion systems since the conventional 
Kondo effect is the formation of a three-body bound state 
between a spin flip and a conduction electron. How¬ 
ever here the three-body wavefunction carries no quan¬ 
tum number and thus is not an order parameter; this is 
why conventional Kondo behavior is associated with a 
crossover and not a true phase transition. 

In the mean-field formulation of hastatic orderj^ a 
spin-1/2 order parameter develops as a consequence of 
a factorization of a Hubbard operator that connect the 
Kramers and non-Kramers states; it is a tensor operator 
that corresponds to the three-body combination 

X acj {R) = |/ 2 a)(/V| = AtW,2,3)V4(l)d(2)V>c(3), 

(159) 

where we have used the short-hand notation 1 = R\ etc. 
and 


A^(ii;l,2,3) = (R 1 ,a;R 2 ,b\X aa (R)\R 3 ,c) 


(160) 


Hastatic order carries half-integer spin and cannot de¬ 
velop via this mechanism. We are then led to the ques¬ 
tion of whether it is possible for Landau order param¬ 
eters to transform under half-integer representations of 
the spin rotation group? At first sight this impossible for 


defines the overlap between the Hubbard operators and 
the bare electron states. In a simple model, this three 
body wavefunction is local, A ^(R; 1, 2, 3) = A ^S(R — 
1 )S(R — 2 )S(R — 3). The factorization of the Hubbard 
operator into a spin-1 fermion and a spin-1/2 boson 


X aa (R)^xl(R)(V*(R)), (161) 

then represents a “fractionalization” of the three body operator. Written in terms of the microscopic electron fields, 
this becomes 

a t c (R\ 1,2,3) 4(l)d'(2)^c(3) = xl(R) (**(R)) • 


This expression can be inverted to give the three body contraction 

1 - 1 - 1 / \ 

V’l(l)V’b(2)V’c(3) = J2 G abc(h2,3-,R)xi(R) (^a(R)) , 

R 


( 162 ) 
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where G^ c (l, 2,3; R) = [Af a c (R; 1, 2,3)]* (see Fig. 14 b) 


(163) 


The asymmetric decomposition of a three-body 
fermion state into a binary combination of boson and 
fermion is a fractionalization process; if the boson carries 
a quantum number, when it condenses we have the phe¬ 
nomenon of “order parameter fractionalization”. Frac¬ 
tionalization is well-established for excitations of low di¬ 
mensional systems, such the one-dimensional Heisenberg 
spin chain and the fractional quantum Hall effeclP^^, 
but order parameter fractionalization is a new concept. 
Unlike pair or exciton condensation, the hastatic order 
parameter transforms under a double-group representa¬ 
tion of the underlying symmetry group, and thus repre¬ 
sents a fundamentally new class of broken symmetries. 
We are currently investigating order parameter fraction- 
alization beyond the realm of URv^SU . The proposed 
three-body bound-state has a nonlocal order parameter, 
and it may be possible to identify a dual theory with a 
local order parameter that breaks a global symmetry. 

B. Experimental Constraints and More Tests 

Let us now return to the situation in URa^SU . As 
we discussed earlier, hastatic order leads to a predica¬ 
tion of a basal-plane moment of order where Tk 
and D are the Kondo temperature and the band-width 
respectively. The transverse moment in our mean-field 
treatment has contributions from both conduction and f 
electrons, and the ratio ^ is very sensitive to the de¬ 
gree of mixed valence of the U ion. Our original cal¬ 
culation assumed 20% 5/ 3 , leading to a predicted basal- 
plane moment of 0.01 \ib- Recent high-resolution neutron 
experiments 72 with momentum transfer along the c- 
axis designed to detect this predicted transverse moment 
have placed a bound p± < 0.0011 pb on the ordered 
transverse moment of the uranium ions, constraining it 
to be at best an order of magnitude smaller than what 
we predicted. 

Clearly we need to reconsider our calculation of the 
transverse moment and understand why it is so small if 
not absent and we are currently exploring a number of 
possibilities: 

• Fluctuations. Amplitude fluctuations of the 
hastatic order parameter are needed to describe 
the incoherent Fermi liquid observed to develop 
at temperatures well above Tho in optical, tun¬ 
neling and thermodynamic measurements IM H 4Q I 56 I 
and they will reduce the transverse moment. We 
note that various pro bes, incl uding X-rays, p- spin 
resonance and NMRp^ 67 75 76 have consistently de¬ 
tected basal plane fields of order 0.5G, consistent 
with the presence of a tiny in-plane moment. 


• Uranium Valence. The predicted transverse mo¬ 
ment is very sensitive to the 5f valence, decreas¬ 
ing with increasing proximity to pure 5/ 2 . More 
specifically it is proportional to the change in va¬ 
lence between Tho and the measurement temper¬ 
ature and thus is significantly smaller than the 
high-temperature mixed valency. It would be very 
helpful to have low temperature probes of the 5f- 
valence. 

• Domains. X-rayf^ muonj^ torque 

magnetometr}!^, c yclot ron resonancd^ and 

NMR measurementthat have indicated 
either a static moment or broken tetragonal 
symmetry were performed on small samples. By 
contrast, the neutron measurements that show no 
measurable moment use large samples 72 74 . The 
apparent inconsistency between these two sets of 
measurements may be due to domain formation 
of hidden order. Such domain structure could 
result from random pinning 78 of the transverse 
moment by defects of random strain fields. The 
situation in UR^SU is somewhat analogous to 
that in S^RuCU, where there is evidence for 
broken time-reversal symmetry breaking with a 
measured Kerr effect and /aSR to support chiral 
p-wave superconductivity, but no surface currents 
have yet been observed. 79 Domains are an issue in 
this system too. 

• x-y order and spin superflow. The current mean- 
field theory has the transverse hastatic vector 

pointing in one of four possible directions at 
each site, corresponding to a four-state clock model. 
The tunneling barrier between these configurations 
is very small. When we expand the effective action 
as a function of </>, the leading order anisotropy will 
have the form 

A E((j>) = U 4 cos40. (164) 

where determines the magnitude of the tun¬ 
neling barrier. Now the anisotropic terms have the 
form e ±24< ^, and since the dependence in T enters 
as e ±l 2 5 the leading dependence of this term on T 
has the form 

E a ~T k |^| 8 (165) 

Now since 'h| 2 ~ y, this implies that the tunnel- 
ing barrier has magnitude 



( 166 ) 
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In our theory we have estimated Tk/D ~ 0.01, so 
that the tunneling barrier is of order 10 -8 times 
smaller than the Kondo temperature. In practice, 
the XY-like basal plane hastatic moments will be 
extremely weakly pinned, with large domain walls 
between Z4 domains, with widths ~ ^ ~ 100 lat¬ 
tice spacings. To our knowledge, such nearly per¬ 
fect XY order, which can lead to spin superflow®^ 
is completely unknown in magnetism: its only 
counterpart occurring in neutral superfluids. This 
opens the interesting possibility that the presence 
of persistent spin currents in the hastatic phase 
jspin oc V0 give rise to a destruction of the stag¬ 
gered moment associated with hastatic order. By 
contrast, near the surfaces, where the tetragonal 
symmetry is broken, the Z 2 pinning is expected to 
be much greater. This might account for why large 
moments and broken tetragonal symmetry only ap¬ 
pear in tiny crystals. 

The central tenet of the hastatic proposal is that Ising 
quasiparticles are associated with the development of hid¬ 
den order and there remain several tests of this aspect 
that can be made. In particular: 

1. Giant Anisotropy in A \3 oc cos 4 0. In this mea¬ 
surement the temperature-dependence of the Ising 
anisotropy of the conduction fluid can be probed to 
confirm that it is associated with the development 
of hidden order. 

2. dHvA on all the heavy Fermi surface pockets. 
Based on the upper-critical field results, we expect 
that the heavy quasiparticles in the a f3 and 7 orbits 
will exhibit the multiple spin zeros of Ising quasi¬ 
particles but to date only the a orbits have been 
measured as a function of field orientation. 

3. Spin zeros in the AFM phase? (Finite pressure) If 
the antiferromagnetic phase is also hastatic, then 
we expect the spin zeros to persist at finite pres¬ 
sures. 


C. Future Challenges 

The obs ervat ion of Ising quasiparticles in the hidden 
order state 22 23 83 84 represents a major challenge to our 


understanding of URu2Si2 ; to our knowledge this is 
the only example of such anisotropic mobile electrons, 
and as we have emphasized, completely unexpected for 
f-electrons in a tetragonal environment. As we have em¬ 
phasized throughout this paper, Ising quasiparticles are 
the central motivation for the hastatic proposal, and a 
key question is whether this phenomenon can be de¬ 
scribed by other HO theories? In particular: 

• Can band theory account for the g(0) observed in 
URu 2 Si 2 ? Recent advances in the understanding of 
orbital magnetizatiorP^^ suggest it may be possi¬ 
ble to compute the g-factor associated with conven¬ 
tional Bloch waves; in a strongly spin-orbit coupled 
system, the orbital contributions to the total energy 
in a magnetic field are significant. It would be par¬ 
ticularly interesting to compare the g{0) computed 
in a density functional treatment of URu2Si2 with 
that observed experimentally. 

• Can other 5 f 1 2 3 theories account for the multiple spin 
zeroes and the upper bound A < 1 K on the spin 
degeneracy of the heavy fermion bands? In partic¬ 
ular, is it possible to account for the observed spin 
zeros without invoking a non-Kramers 5 f 2 doublet? 
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